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O : ABSTRACT 

o 

■ We present results from RXTE data obtained during a systematic moni- 

[Jh \ toring program of four Seyfert galaxies (NGC 5548, NGC 5506, MCG -6-30-15 

' and NGC 4051). We studied the variability of three hardness ratios derived 

^ ■ from the light curves in four energy bands {HRl, which describes the continuum 

variations, and HR2, HR3 which are sensitive to the iron line and reflection com- 
ponent variations with respect to the continuum, respectively). All the objects 
Oj^ ■ show similar spectral variations in all ratios. In order to interpret the results we 

^ ■ computed the hardness ratios corresponding to a simple spectral model of a power 

^ ' law plus iron line plus reflection component. In order to derive the model HR2 

O ■ and HR3 colors, we considered two possibilities: a) variations with constant line 

^ ■ equivalent width and reflection parameter R (the case of a reflecting/reprocessing 

material that responds with a short delay to the continuum variability) and b) 
variations with constant line and reflection flux (the case of a reprocessor that 
c/3 . does not respond to the fast, intrinsic variations). The overall, mean observed 

trends can be explained by spectral slope variations (AF ~ 0.2 — 0.3, and ~ 1 
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for NGC 4051), and a constant flux Fe line and reflection component, although 
the existence of a line component which is variable on short time scales cannot 
be excluded. Finally, we flnd that the data are not consistent with an increase 
of R with flux for individual sources, indicating that, as a single source varies, 
softer spectra do not correspond to larger R values. 
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Subject headings: Galaxies: active - Galaxies: nuclei - X-rays: Galaxies - Galax- 
ies: individual(NGC 4051, MCG -6-30-15, NGC 5506, NGC 5548) 

1. Introduction 

X-ray emission is a ubiquitous property of Active Galactic Nuclei (AGN). The spectral 
characteristics of this emission have been studied in detail the past decade and our under- 
standing of the X-ray emission properties has improved considerably. In general, the X-ray 
spectra of AGN have a power-law form over a large range in energy. The photon spectral 
index has a typical value of F ~ 2 (Nandra & Pounds 1994). Few sources also show a high 
energy cutoff which can be modeled as an exponential cutoff with e-folding energy of a few 
100 keV (Gondek et al. 1996). Furthermore, the presence of the 6.4 keV Fe line and the 
so-called "Compton reflection hump" indicates the existence of cold, optically thick material 
that reprocesses and reflects the medium energy X-rays. It is generally believed that this 
material is located near the primary X-ray source (i.e. the accretion disc). However, this 
is not the only possibihty. The cold material could also be in the form of a torus which is 
surrounding and obscuring (depending on the inclination) the central nucleus. The existence 
of this torus is predicted by models that unify Seyfert 1 and 2 galaxies (e.g. Antonucci 1993). 
The contribution of the torus to the spectrum of Seyfert 1 galaxies (by scattering part of the 
primary radiation into the line of sight) can in principle explain the reflection bump observed 
in the spectra of AGN (Ghisellini, Haardt & Matt 1994). 

Based on the detailed study of the X-ray spectra of the Seyfert 1 nuclei the currently- 
popular scenario (the "standard picture" say) is that X-rays from AGN are produced by 
Compton upscattering of soft photons in a hot "corona" located above the accretion disc (e.g. 
Haardt & Maraschi 1991, 1993). The X-ray continuum illuminates the disk, which is assumed 
to be neutral, producing optical/UV emission via reprocessing, the iron line and the reflection 
hump. Apart from the energy spectra, the X-ray flux and spectral temporal variations can 
also constrain the emission mechanism. In fact, the spectral variations observed during long, 
high signal-to-noise observations of NGC 5548 and NGC 7469 (Petrucci et al. 2000, Nandra 
et al. 2000) have provided strong support to the Comptonization hypothesis as the origin 
of X-rays in AGN. In the former case, during a long Beppo SAX observation the spectrum 
became "harder" (i.e. flatter) and the corona temperature decreased together with the 
[2 — 10] keV flux whereas the total luminosity remained essentially constant. In NGC 7469, 
simultaneous RXTE and lUE observations revealed strong X-ray spectral changes which 
were correlated with the UV flux. In both cases, the observed spectral variations can be 
naturally explained in terms of thermal Comptonization models since an increase in the UV 
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photons reduces the temperature of the hot corona and produces a "softer" (i.e. steeper) 
X-ray power law slope. 

In this work we present the results from a spectral variability analysis of four Seyfert 
galaxies, namely NGC 4051, MCG -6-30-15, NGC 5506 and NGC 5548. These objects 
have been monitored with RXTE on a regular basis since 1996. The main purpose of the 
observations was to compute accurately the shape of the power density spectrum over a 
broad range of time scales in order to reveal the existence of characteristic time scales (low 
or high frequency "breaks"). These time scales can in principle provide model- independent 
black hole mass estimates through their comparison with similar time scales that have long 
been identified in the power spectra of Galactic X-ray binaries. The results of this study are 
presented elsewhere (Uttley, M'^Hardy & Papadakis 2001). Here we focus on the variability 
properties of the energy spectra of the sources. 

All the sources are strongly variable in the X-rays. MCG -6-30-15 and NGC 5548 are 
classified as Seyfert 1 nuclei, NGC 5506 is a Seyfert 2 galaxy while NGC 4051 is a Narrow 
Line Seyfert 1 (NLSl) galaxy. They are typical examples of their class and taken as a whole 
they can be considered as a representative group of the nearby Seyfert galaxies. 

There are two main advantages of the RXTE light curves that we study in the present 
work: i) they are of high quality (i.e. high signal to noise ratio) and ii) they span regularly 
a long time interval (~ 3 years for each source). As a result, we are able to study spectral 
variations covering the highest, lowest and intermediate fiux states over a few years. This 
is not possible to achieve with the typical X-ray observations of AGN which last for ~ few 
days (maximum). 

Using the same RXTE data. Lamer, M'^Hardy & Uttley (2000) and Lamer et al. (2001) 
have already presented a study of the spectral variability of NGC 5506 and NGC 4051. Their 
study was based on model fitting of flux-averaged energy spectra. In this paper, we follow a 
different approach. Typically, each source was observed for ~ 1 ksec during each pointing. 
Although this is not long enough for a direct study of the energy spectrum it allowed us 
an accurate fiux measurement in the [3 — 15] keV band. Using the RXTE light curves at 
different energy bands we computed hardness ratios, produced color-flux diagrams and used 
these diagrams to investigate the properties of spectral variations. As a result we are able to 
study the spectral variations in a systematic way, i.e. at each individual observation without 
considering broad "flux states" for each source. In this way, we cannot investigate whether 
a given model describes well the energy spectrum of the sources. However, using the simple 
"power law -|- iron line -|- reflection component" model which provides a good representation 
of the X-ray spectrum of AGN, we were able to study in detail the relationship between the 
continuum, line flux and reflection component. 
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The paper is organized as follows. In the following section we present our data and in 
Section 3 we present the results from the analysis of the variability of the hardness ratios. In 
Section 4 we use the standard "power law plus reflection component plus narrow iron line" 
model to compare quantitatively the observed trends in the color-flux diagrams with specific 
hypothesis on the variability properties. In Section 5 we compare our results with previous 
studies, and in Sections 6 and 7 we present the discussion of our results and our conclusions. 

2. Observations and Data reduction 

We use data obtained with the PCA on board RXTE covering a three year period during 
observing cycles 1-3 when the PCA gain setting was constant. Details of the data sampling 
scheme and of the time of the first and last observation of the sources can be found in Uttley 
et al. (2002). Briefly, the four sources were observed ~ 120 times during the period from 
April/May 1996 until December 1998/February 1999. Each observation was ~ 1 ksec long 
and the observing scheme was designed so that the broadest range of time scales would be 
covered with the minimum observing time. The objects were observed once every week in 
the first year of the monitoring campaign and then every two weeks. During some periods 
the monitoring was more intensive (twice daily or daily for a period of two up to four weeks) 
to sample the shorter time scale variations. The data were reduced using FTOOLS V4.2. 
We used data from the top layer of the PCUs 0,1,2 only and estimated the background light 
curves using the L7 background model (for details see Uttley ct al. 2002). 

The background subtracted [2 — 10] keV light curves of the four sources, normalized 
to their mean, are shown in Figure 1. All are significantly variable. The peak to peak 
variability amplitude is ~ 40,5,3 and 4 for NGC 4051, MCG -6-30-15, NGC 5506 and 
NGC 5548 respectively. Therefore NGC 4051 stands out as the most variable source as is 
typical of NLSl nuclei (e.g. Leighly 1999a). During the iiXTE observations, NGC 4051 went 
through an unusually "low-state" when the source reached a very low flux level (the period 
between Day 600 and 750 in Figure 1). In fact, there is the possibihty that the central source 
even "switched-ofF' , leaving only the reflection spectrum from cold matter to be detected 
(Guainazzi et al. 1998, Uttley et al. 1999). For that reason, we excluded the points in the 
NGC 4051 hght curves between Day 600 and 750 from the data analysis. 

Finally, the RXTE observations of NGC 5548 are contaminated by the nearby BL Lac 
object IE 1415.6-1-2557. Chiang et al. (2000) have estimated that the contaminating contri- 
bution of this object to the 2 — 10 keV NGC 5548 light curves (using the same three PCUs 
as ours) is ~ 2 counts/sec with a modulation smaller than 0.8 counts/sec. Since we observe 
a peak to peak count rate difference of ~ 20 counts/sec, the IE 1415.6-1-2557 contribution 
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to the overall variability properties of NGC 5548 should not be crucial. 

3. Data Analysis 

3.1. Hcirdness Ratios 

We extracted hght curves in the following energy bands: [3 — 5], [5 — 7], [7 — 10] 
and [10 — 15] keV. These intervals were chosen so as to separate as much as possible the 
various emission components in the spectrum of the objects. The first and third band should 
be representative of the primary continuum mainly, while the iron line and the refiection 
component emission should contribute in the [5 — 7] keV and [10 — 15] keV bands respectively. 
From these light curves we calculated three hardness ratios: HRl = Ci^iQkeV /Cz-bkeVi 
HR2 = Cj^iokev/Cb-Jkev and HR3 = Cio-isfcey/Cr-iofceV iCEi-E2 represents the count rate 
in the energy band El — E2). These are sensitive to variations : i) of the continuum shape 
(HRl), ii) of the continuum to the iron line ratio {HR2) and iii) of the reflection component 
to the continuum ratio {HR3). 

Figure 2 shows a plot of the hardness ratios as a function of time since the beginning 
of the RXTE observations for each source. The sohd fines show the HR weighted average 
values {HR). These values, together with the root mean square variability amplitude (i.e. 
(Jrms) of tfie HR light ciuvcs are reported in Table 1. A constancy test shows that all the 
HR values exhibit statistically significant variations, except from the HR?) plot of NGC 5548 
for which the significance is only marginal. The HRl light curve of NGC 4051 has the largest 
arms among the HRl light curves. This source shows the largest amplitude flux and spectral 
variations. 

The main conclusion from Table 1 is that the mean HR values are roughly similar 
between the different objects. This implies that the mean energy spectrum has approximately 
the same spectral shape in all sources. However, there are also some differences. For example, 
NGC 4051 shows the smallest HR values, indicating that its energy spectrum is softer than 
the spectrum of the other sources. This is consistent with results from previous studies which 
have showed that the X-ray energy spectra of NLSls are "softer" than the spectra of typical 
SI galaxies (e.g. Leighly 1999b). The HR values of NGC 5506 (the only S2 source in our 
sample) are similar to the respective values of NGC 5548 and MCC -6-3015 (the SI sources 
in the sample) except from HRl. Its large value is probably due to the presence of a large 
absorbing column towards that source which causes the energy spectrum to appear "harder" 
than the spectrum of other sources. 
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3.2. Color-Flux diagrams 

Having established the existence of significant spectral variability in all sources, we 
then investigated whether the spectral variations are correlated with the source flux. Figure 

3 shows the plot of HRl as a function of the normalized ([3 — 5] + [7 — 10]) keV count 
rate {C(^3_5keV+7-iokeV)„orm.)- Figures 4 and 5 show a plot of HR2 and HR3 as a function 
of the normalized [3 — 5] keV count rate {C(3-5kev)„orm.)- 1^ Figure 3, we have chosen the 
([3 — 5] + [7 — 10]) keV count rate as representative of the continuum flux, since the sum 
of the source signal in the two bands should minimize the possibility to introduce artiflcial 
correlations due to the inter-dependency of the HRl and C(3_5+7_iofcey) variables. In Figures 

4 and 5, we have used only the [3 — 5] keV count rate since both HR2 and HR3 are 
independent of C(3_5fcet))- In all cases, we have normalized the count rates to their mean in 
order to compare the color-flux plots of the four sources. 

The points in the color-flux plots of all the sources cluster in a well deflned, "continuous" 
region rather than forming a scatter diagram or fllling separate "islands". Although it is 
not clear from the power density spectra of the sources whether we have observed their 
maximum/minimum flux states with the present iiXTi? observations (see discussion in Uttley 
et al. 2002) we are certain that we have observed all the "intermediate" flux states of each 
source between their largest and lowest flux state in the RXTE light curves. 

The color-flux trends are similar for each object: HRl systematically decreases (i.e. the 
energy spectrum softens) as the flux increases. A similar trend is observed for HR?,. On the 
other hand, the HR2 variations appear to be independent of the source flux. 

In order to quantify the above relationships we fitted the data for each galaxy using 
a power law function of the form: HR = aiC^^^. We also tried a linear relationship (i.e. 
HR — ai + hiX Cnorm) but the goodness of fit was worse in all diagrams. The model fitting 
was done taking account of the errors in both variables. The best fitting results are listed in 
Table 2. The best fitting hi values in the HR2/C{z-bkeV)^^^^ diagrams are consistent with 
zero (contrary to the HRl/C(^3^5keV+7-iokeV)^^^^ s-nd HR'i/C{z-5keV)^^^^ diagrams). This 
result shows that the HR2 variations are indeed independent of the source fiux. The only 
exception is NGC 4051 which shows a HR2/C{^^^heV)^^^^ slope which, although small, is 
significantly different from zero. 

NGC 4051 also shows the largest absolute hi values in all diagrams. The spectral 
variations in this source are thus a "steeper" function of the fiux than in the other sources, 
i.e. for the same amplitude flux variations, the spectral variations in NGC 4051 are larger. 
Furthermore, the NGC 5506 hi value in the H Rl / C (^-^keV+i-wkeV)^^^,^ diagram is different 
from the values of NGC 5548 and MCG -6-30-15. This difference cannot be attributed to the 



- 7- 



presence of a large absorbing column towards NGC 5506, unless the column density varies 
with the X-ray source flux. However such behavior has never been observed in S2 nuclei 
(Risahti, Elvis and Nicastro 2001). 

As Figures 4,5 and 6 show, a power law describes well the overall trend in all the 
color-flux plots. However, it does not give give a statistically acceptable fit, except for the 
H R3/ C(3^5keV)norm dlagrams (this is probably due to the larger errors associated with the 
HR3 values). There is significant scatter around the mean trends which shows that the 
spectra are significantly variable at any given "fiux state" of the sources as well. 

The spectral variability behavior of the sources could be studied with the frequently 
used color-color diagrams as well. For example, HR2 vs HRl diagrams could be used to 
investigate, directly, whether the variations of the line's EW are correlated with the primary 
spectrums slope changes. However, as we show in the Appendix, the fact that both HRl and 
HR2 use the count rate in the [7 — 10] keV band, implies that there could exist misleading 
correlations in the HR2/HR1 diagram. For that reason, we do not use any color-color 
diagrams to investigate the spectral variability behavior of the sources. 

4. Determination of spectral parameters from the hardness ratios 

In order to understand the spectral variability behavior of the four sources in terms of 
commonly used spectral parameters, we used the widely used model consisting of a power 
law plus Fe line plus a reflection component. This model flts well the average X-ray/7-ray 
spectrum of Seyfert 1 nuclei, especially in the 2-20 keV energy range (Nandra & Pounds 
1994). We computed the model hardness ratios under different assumptions and compared 
them with the data as described in detail below. This comparison allows us to derive the 
main spectral characteristics (i.e. spectral index, reflection and iron line equivalent width) 
and to discuss the spectral variability exhibited by the different objects in our sample. 

4.1. The method 

The model consists of a power law continuum of the form: /e = NE~^, where N is the 
normalization at 1 keV. We included a Compton reflection component using the PEXRAV 
model (Magdziarz & Zdziarski 1995) in XSPEC. This extra component calculates the X-ray 
spectrum when a source of X-rays is incident on optically thick, neutral (except hydrogen 
and helium) material. We fixed the inclination angle for the disk at i = 30° (the shape of the 
refiection spectrum below 20 keV is relatively independent of the inclination angle). The iron 
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and light element abundances were kept fixed at the solar abundance values. The strength 
of the reflection component is governed by the parameter R, representing the strength of the 
reflected signal relative to the level of the incident power-law continuum. We also added a 
narrow Gaussian line {(June = 0.1 keV) to simulate the iron line. The line energy was kept 
fixed at 6.4 keV so that the line is completely characterized by its equivalent width {EW). 
Finally, except from NGC 5506, we do not consider the effects from neutral absorption in our 
calculations, since the absorbing column towards the direction of the sources is consistent 
with the Galactic value and does not affect the X-rays at energies > 2 keV. For NGC 5506 
we fix the value of Nh at 3.6 x 10^^ cm^^ (Perola 1998). 

Assuming that this model gives a good representation of the spectral shape of the 
sources, we have a one-to-one correspondence between the data (the count rate in the four 
energy bands) and the model parameters {N, F, EW and R). As a result, the observed 
count rates can be used to calculate the model parameters as described below. 

Using XSPEC (vll.O), we computed the model count rate expected in different en- 
ergy bands from the power law continuum, the refiection component and the Iron line for 
different values of the spectral index F, but fixing N — 1, R — 1, and EW — 250 eV 
(Csi-E2(r), C|i_s2(F) and C|rs2(F), respectively, where {EI-E2) = [3-5], [5-7], [7-10] 
and [10 — 15] keV). Then, for each observation i, the observed count rate at each energy 
band is given by. 



Ce\-e2 — Ni 



(1) 



We have thus a system of four equations with four unknown variables (A^j, Fj, Ri and 
EWi) which is easily solvable for each pointing. 

Since we are interested in the overall spectral behavior of the four sources during the 
total observational period, we did not solve the system of equations (1) for each observation. 
Instead, we used these equations to calculate the spectral slope and normalization at various 
flux states along the best fitting lines shown in Figure 3. Then we computed the HR2 and 
HR'i model ratios considering two different possibilities for the response of the reflecting 
material to the continuum variations: (a) The material responds with a short delay to the 
continuum variability, in which case the line's EW and the reflection component's parameter 
R should remain constant during the spectral variations. This is expected in the case when 
the material is located close to the X-ray source (i.e. the accretion disk), (b) The reflect- 
ing/reprocessing material does not respond simultaneously to the fast, intrinsic variations, 
as expected if the reprocessor is located away from the central source (i.e. the obscuring 
torus). In this case the line and the refiection component flux should remain constant. 
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The model HR2/C^s^5keV)„orm HRS/C^s^^keVUorm curves, in the case of constant 
EW and R, are plotted in Figures 4 and 5. In Figure 4, the dashed, dotted and dashed/dotted 
lines show the model curves in the case of constant EW = 0, 250 and 500 eV respectively 
{R is kept constant to 1). The dashed, dotted and dashed/dotted lines in Figure 5 show the 
model curves in the case of constant R = 2,1 and respectively. 

To produce the model diagrams in the case of a remote reflector, we used the average 
HR values, listed in Table 1, to find the spectral parameters (i.e., F, A^, EW and R) of 
the time average continuum of the four sources. Using XSPEC, we then deduced the flux 
of the line and of the reflection component in the average spectra. These fluxes were kept 
constant in order to estimate the model HR2 and HR3 values, at each C(3_5fcey)„orm) using 
the respective (F, N) values (computed for R — 0). The corresponding model curves are also 
plotted in Figures 4 and 5 (solid hues). 

4.2. Results 

4-2.1. The average spectra 

Using equation (1) and the mean H~R1, HR2 and HR?, values of each object, we esti- 
mated the mean spectral index (F), mean reflection normalization (i?) and mean iron line 
equivalent width [EW). The results are listed in Table 3. The values in this table show 
that NGC 4051 possesses the largest average R and EW values among the four sources. The 
mean spectral characteristics of the Seyfert 2 NGC 5506 are similar to those found for the 
two SI objects, NGC 5548 and MCG -6-30-15. 

Our results are in good agreement with the respective values that have been reported in 
the past, based on model fltting of the full band energy spectra. For example. Lamer et al. 
(2000, 2001), flnd an average value of F ~ 2.1, R ~ 1.2 and EW ~ 300 eV for NGC 5506, 
and F ~ 2.05, EW ~ 400 eV for NGC 4051. These values are consistent with the values 
listed in Table 3 for the same objects. The average values for the hne's EW that have been 
reported in the past for MCG -6-30-15 and NGC 5548 (Wang et al. 1999, Lee et al. 1999, 
and Chiang et al. 2000) are in agreement with the respective values in Table 3. The average 
it! value that we find for MCG -6-30-15 also agrees with the values that Lee et al (1999) 
report, while it is slightly larger in the case on NGC 5548 (0.8 as opposed to ~ 0.5 in Chiang 
et al 2000). Finally, the slope values that we flnd for MCG -6-30-15 and NGC 5548 are in 
good agreement with the values reported in the past for the same objects. 

Interestingly, we find a correlation between (i?) and (F) as well as a correlation between 
{EW) and (F). When comparing the average spectra of the four sources, we see that the 
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sources which show harder spectra show smaller (average) R and EW values as well. This 
correlation is in agreement with the claims of Zdziarski et al. (1999). 

4.2.2. The primary continuum variability 

The dashed lines in Figure 3 show the best fitting curves to the observed color-flux 
diagram of the four sources. On these hues, we have plotted the spectral slope values for a 
few normalized count rate values. These slope values are computed as explained in Section 
4.1 in the case of i? = 1. According to the model predictions, the observed HRl variations 
correspond to significant continuum slope changes. Differences in the spectral index up to 
a Ar ~ 0.3, ~ 0.2 and ~ 0.2 are necessary to explain the HRl variability in MCG -6-30- 
15, NGC 5506 and NGC 5548 respectively. In NGC 4051, differences larger than AF ~ 1 
are needed to explain the observed HRl variations. Although we have excluded from our 
analysis the points that correspond to the "off-state" of this source, Uttley et al. (1999) have 
shown that the decline of the source towards this "off-state" is gradual. If the spectrum in 
this state is simply a reflection component produced by the torus and remains present at all 
times, then this is probably the dominant component in the X-ray spectrum of NGC 4051 
in the observations that have the lowest count rate in Figure 3. As a result, the spectral 
slope that we measure could be affected by the shape of the reflection component, and the 
primary slope could be steeper than what we estimate. In fact. Lamer et al. (2001), find 
slope variations between F ~ 1.6 and 2.3 when they subtract from the source's spectra at 
the various flux states the "off-state" spectrum. As expected, their upper limit is consistent 
with ours (~ 2.4, see Figure 3) while their lower limit is larger than ours. 

In order to investigate further the effects of a constant ffux reffection component on 
the spectral slope variability, we used the flux of the reflection component in the average 
spectrum of the sources (see Section 4.1) and assumed that it remains constant. In this case, 
since the reflection component contributes to the [7 — 10] keV band more than it does to the 
[3 — 5] keV band, it can result in an artificial fiattening of the photon index, mainly when 
the source is at low flux states. For that reason, taking into account the constant reflection 
component flux at each energy band, we computed again the [A^, F] values (using the best 
fltting line to the HRl vs flux plot, as explained in Section 4.1). Figure 6 shows again the 
plot of HRl as a function of the normalized ([3 — 5] + [7 — 10]) keV count rate, with the 
HRl values of NGC 4051 during the "low-state" included as well (in order to examine the 
intrinsic spectral slope values when the source is dominated by the refiection component). 
On the best fitting lines, we mark the new spectral slope values for a few normalized count 
rates, like we did in Figure 3. 
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As expected, due to the presence of a constant flux reflection component, smaller AF 
variations are necessary in order to account for the observed HRl variations, and the "bound- 
ary" r values are now steeper. However, the differences between the results shown in Figures 
3 and 6 arc small. The assumption of a constant flux reflection component alone cannot ac- 
count for the observed HRl variations. Intrinsic spectral slope variations of the order of 
~ 0.2 are still needed in order to explain the spectral variations in MCG -6-30-15, NGC 5506 
and NGC 5548. Even in NGC 4051, intrinsic variations of F between ~ 1.5 and ~ 2.4 are 
still implied by the HRl variations. In fact, despite the presence of the constant flux reflec- 
tion component, a primary component (albeit with a small normalization) with a very flat, 
variable slope needs to be present in order to account for the HRl values when the source 
is in the low/off-state. 

More important effects on the implied F variations can have the variability behavior of 
the broad red wing of the Fe line in MCG -6-30-15. Based on ASCA observations, Iwasawa 
et al. (1996) found that the iron line in this source consists of a narrow core at an energy 
of ~ 6.4 kcV, and a broad wing extending to below 5 keV. While the narrow component 
cannot contribute significantly to the count rate in either the [3 — 5] or [7 — 10] keV band, the 
broad red wing will contribute to the observed [3 — 5] keV count rate. Iwasawa et al. (1996) 
found that the broad component can be fitted well with a Gaussian centered at 5.5 keV with 
a dispersion of 0.64 keV. They also found that the EW of this component component is 
variable, decreasing from ~ 600 eV to ~ 100 eV with increasing flux. Consequently, when 
the source is at low state, the broad wing can contribute significantly to the [3 — 5] keV band 
count rate, resulting in a flatter intrinsic slope. 

In order to investigate this possibihty, we recalculated the [iV, F] values for MCG -6- 
30-15, assuming that the EW of the fine's red wing decreases from ~ 0.6 keV to ~ 0.1 
keV, when the source is at the lowest /highest states in the present RXTE light curve. As 
expected, when the source is above its average flux state, the presence of the broad iron line 
wing does not affect significantly the estimation of the spectral slope value. However, as the 
flux decreases, the slope values that we estimate become flatter. At the lowest flux state (i.e. 
when C(3_5fcey+7-iofcey)„„^^ ~ 0.4 in the MCG -6-30-15 plot of Figures 3, 6) we flnd a value 
of F ~ 1.7, as compared to F ~ 1.9 that we flnd when we do not consider the effects of the 
fine's red wing. We conclude that, if the variabifity behavior of the extended red wing of 
the Fe line in MCG -6-30-15 on long time scales is similar to the behavior exhibited during 
the ASCA observations, the observed HRl variations of this source imply variations of the 
primary spectral slope of the order of AF ~ 0.4 — 0.5. 

Finally, using the [iV, F] values at each flux state as shown in Figure 3, we computed 
the primary power law flux in the energy range between 1 — 300 keV. Figure 7 shows a plot 
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of this flux as function of the normahzed source count rate. In NGC 4051, we find that 
the luminosity variations have an amplitude much smaller than that of the observed flux 
variations. In this source, the large amplitude [3 — 15] keV flux variations appear to be 
caused mainly by the large amplitude continuum slope variations. In the other sources, the 
peak to peak amplitude of the luminosity change is ~ 2. This is smaller but comparable with 
the amplitude of the observed flux variations. However, since the variability behavior of the 
high energy cut-off component cannot be studied with the present RXTE data, it is possible 
that the amplitude of the luminosity variations is smaller in reality if there are significant cut 
off variations. In fact, Petrucci et al. (2000) find that in NGC 5548 (the source that shows 
the largest amplitude luminosity variations in Figure 7) the cut-off variations are associated 
with the spectral slope variations in such a way so that the total X-ray luminosity remained 
constant during the long Beppo SAX observation. 



4-2.3. The line variability 

Filled squares in Figure 4 show the average HR2 values over various flux (i.e. C(^3-5keV)norm) 
bins. The size of the bins was variable so that 15 points were included in each one. These 
points indicate the average HR2 behavior as a function of the soft band flux and are close 
to the best fitting lines of the HR2/C(^3^5keV)^^^^ diagrams. 

The "constant EW" model predictions do not agree well with the average HR2 values. 
Instead, the model curves suggest rather complex variations for the line's EW in all sources. 
It appears that, on average, the line's EW decreases as the source flux increases. For 
example, for NGC 5548, when the source is at the lowest flux state the line's EW is ~ 250 
eV. As the source flux increases, the line's EW decreases to ~ 50 eV. In NGC 4051, when the 
normalized soft band count rate is below ~ 0.7, the line's EW appears to increase strongly, 
reaching a value larger than 500 eV at the very low flux states. The tendency of the line's 
EW to decrease with increasing source flux also appears in the color-flux diagrams of NGC 
5506 and MCG -6-30-15 but the amplitude of the EW variations is smaller in these two 
sources. 

In the same Figure, we have also plotted the "constant line flux" model curves. In this 
case, the model predictions are consistent with the average HR2 values, except perhaps NGC 
5506 where the model curve hes above the average HR2 values. We can get an agreement 
between the model predictions and the NGC 5506 data if we increase slightly the flux of the 
line that we add to all the model spectra. 

Column 2 and 3 in Table 4 hsts the slope of the model curves shown in Figure 4. The 
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bcEW values correspond to the slope of the model curves under the assumption of constant 
line's EW, while bcLF is the slope of the constant line flux model curves. Comparison 
of these values with the best fitting slope values listed in Table 2, shows clearly that the 
constant flux model predictions agree better with the data; their slopes are consistent with 
the best fitting slope values. On the other hand, the difference between the constant EW 
model curve slopes and the best fitting slopes is larger than 3a in all sources. 

4.2.4- The reflection component variations 

As before with Figure 4, filled squares in Figure 5 show the average HR3 values as a 
function of the soft band fiux. In the same Figure, we also show both the "constant R" and 
"constant refiection flux" model curves. 

Due to the large scatter in the H R3 / C (j>,^r^kev) diagrams of all sources, and the 
similarity of the model curves under the two different assumptions, it is not clear which 
model curves flt better the data. In the MCG -6-30-15 and NGC 5506 plots, comparison 
between the constant R model curves and the average HR3 values show that R slightly 
decreases with increasing source flux, in agreement with the EW dependency on the flux in 
these sources. However, most of the average HR?, points are close to the constant R — 1 
model curve, and only the lowest flux HR?, value is clearly closer to the constant R — 2 
curve. In NGC 5548, R appears to remain constant during the spectral variations at a value 
smaller than 1. Similarly, in NGC 4051, the average HR3 values imply spectral variations 
under constant R (with a value between 1 and 2). Only the lowest flux points suggest a 
decreasing R value with decreasing source flux. However, this is again caused by the fact 
that the spectral slope that we measure from observations with the lowest flux represents 
the shape of the reflection component itself. Hence the F values are larger and the R values 
are smaller than in reality (see discussion in section 4.2.1) 

The constant reflection flux model curves also provide a good flt to the average HR3 
values in all sources. Columns 4 and 5 in Table 4 list the slopes of the constant R and 
constant reflection flux model curves {bcR and bcRp respectively). Comparison between 
these values and the best fitting slopes of the H R3/ C(^3_5keV)^^^^ diagrams (listed in table 
2) shows that all model curves are consistent with the data. The only exception is the bcR 
value of NGC 5506, which show a 5a difference with the best fitting slope. 
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5. Compcirison with previous studies 

The spectral variability properties of the four sources that we study in the present work 
have already been addressed previously. In most cases, the results are based on spectral 
fitting of the power law plus reflection model to the observed energy spectra. 

All studies have shown that as the source brightens the continuum spectrum steepens. 
Differences in the spectral slope up to AF ~ 0.3, 0.2 and 0.15 have been observed in MCG - 
6-30-15 (Lee et al. 1999), in NGC 5506 (Lamer et al. 2000) and in NGC 5548 (Chiang et 
al. 2000; Petrucci et al. 2000). Our results are consistent with the previous studies. We 
find that the observed HRl variations correspond to larger spectral slope changes. This is 
because we have detected flux variations of a larger amplitude. In the case of NGC 4051 we 
find a larger AF variations compared to Lamer et al. (2001), but as we discussed in section 
4.2.1, this is due to the significant contribution of the "off-state" spectrum to the source's 
spectra at low flux states. 

Variabihty of the iron hne's EW and of the reflection fraction R have also been observed. 
Using ASCA data of 39 AGNs, Nandra et al. (1997) found a clear decrease in the iron line's 
EW with increasing source luminosity. Lee et al. (2000), Lamer et al. (2000) and Chiang 
et al. (2000) also find that the iron line's EW decreases with increasing source flux in 
MCG -6-30-15, NGC 5506 and NGC 5548 respectively. These results are consistent with 
ours. Furthermore, the dependence of the combined "narrow" and "disk" line's EW on the 
source's flux level that Lamer et al. (2001) find in NGC 4051 is entirely consistent with 
our results shown in Figure 4. On the other hand, Wang et al. (1999) have observed a 
positive correlation between the hne's EW and the source flux in NGC 4051 using ASCA 
data, contrary to our results. We note however that NGC 4051 was "weakly" variable (factor 
of 2) during the ASCA pointing in comparison to variations of more than a factor 10 in our 
case. Furthermore our H R2 / C(^3-^keV)^^^^ diagram for this object shows large scatter around 
the mean trend. We thus believe that the results of Wang et al. (1999) are consistent with 
our conclusions. 

Lamer et al. (2000) find that the reflection fraction in NGC 5506 decreases with source 
flux; in particular it is not correlated to the continuum slope as derived from direct fitting 
of RXTE spectra. Chiang et al. (2000) find that the reflection fraction remains constant 
in NGC 5548. These results are also consistent with ours. On the other hand, Lee et al 
(2000) find that the reflection fraction in MCG -6-30-15 increases with source flux (and 
is anticorrelated with the hne's EW) contrary to what we obtain. Finally, in NGC 4051, 
Lamer et al. (2001) find that the reflection fraction of the "disk component" (i.e. the primary 
spectrum without the contribution of the "off-state" spectrum) is very small. Therefore if the 
only reflection component is the one that is observed in the "off-state" then its flux should 
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remain constant. As Figure 5 shows, our results are entirely consistent with this possibility. 

6. Discussion 

We have studied in detail the spectral variability of four Seyfert galaxies (NGC 4051, 
MCG -6-30-15, NGC 5506 and NGC 5548) over a ~ 3 years period. Using RXTE data we 
computed hardness ratios which are sensitive to the primary slope, iron line and reflection 
component variations (HRl, HR2 and HR3 respectively). Our results show that all the 
sources exhibit significant spectral variability. We find that HRl and HR3 decrease as the 
source flux increases while the HR2 variations are independent of the source flux. The 
spectral variability behavior of the Narrow Line Seyfert 1 galaxy NGC 4051 and the Seyfert 
2 galaxy NGC 5506 is similar to the behavior of MCG -6-30-15 and NGC 5548, the two 
Seyfert 1 nuclei in our sample. There are also minor differences. NGC 4051 in particular 
shows the largest spectral slope variability amplitude among the four sources. 

Our approach to study the spectral variability of the sources is different from previous 
studies. Instead of a detailed spectral study that uses the full energy spectrum (usually in 
the ~ 2 — 15 keV band with RXTE data) at a few epochs (at best), we use a large number 
of short exposures over a long period, and compare the observed color-flux diagrams with 
the predictions of the simple "power law plus Compton reflection component plus iron line" 
model. The hardness ratios cannot be used in order to investigate whether the model can fit 
well the X-ray spectrum of the sources or not, however, it is known from past studies that 
this model does provides a good representation of the X-ray energy spectra of AGN. In this 
case, due to the large number of observations which are regularly distributed over a long 
period (much longer than the typical duration of the AGN observations) the hardness ratios 
can be used to investigate, accurately, which components in the X-ray spectra of AGN are 
responsible for the correlations that we find in the color-flux diagrams. 

We did this by constructing various model color-flux curves. For the HR2/C(^s-5kev)„orni 
and HR3/C(^3^^kev)nor,n iiiodel curves we considered two possibilities for the reflecting and 
the line producing material. First, we produced model curves assuming that the material 
responds fast to the continuum variations. In this case the hue's EW and the reflection 
parameter it! remain constant during the spectral slope variations. Alternatively, in the case 
when the material does not respond quickly to the continuum variations, it is the flux of the 
line and the reflection component that remains constant. 

Our main results from the comparison of the observed color-flux plots with the predic- 
tions of the "power law -|- reflection -|- iron hue" model are the following: 
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1) The primary continuum slope is variable. Averaging HRl in fixed flux intervals, the 
derived slope is positively correlated with the soft band flux. At each flux state, there are 
also significant slope variations around the mean slope that corresponds to that state. 

2) As the soft band flux increases, the line's EW decreases (Figure 4). There is an 
indication that the reflection parameter R also decreases with increasing flux in the case of 
NGC 5506 and MCG -6-30-15 (Figure 5). In NGC 4051 and NGC 5548, R remains constant 
during the primary spectral and flux variations. In any case, the HR3 color-flux diagrams 
do not show any increase of R with increasing source flux (i.e. as the primary spectrum 
steepens). 

6.1. Power- law slope variability 

Spectral slope variability is predicted by thermal Comptonization models. An increase 
of the soft input photon flux will result in a reduction of the corona temperature, a steepening 
of the intrinsic spectrum and hence an increase in the soft X-ray band flux. This response of 
the X-ray spectrum to the soft input photon changes has already been observed in NGC 7469 
(Nandra et al. 2000) and NGC 5548 (Petrucci et al. 2000). In the simplest case, we expect 
a correlation between the spectral slope and the soft band flux, as observed in Figure 3. It 
is remarkable how similar the average spectral behavior is for three of the sources and for 
all the intensity states covered. In particular 5506, a Seyfert 2, follows the same behavior as 
5548 and MCG -6-30-15. 

The fact that the spectral slope variations of NGC 4051 have a larger amplitude could 
be due in part to a larger variation of the strong soft component known to be present in this 
source, affecting mainly the [3 — 5] keV band, hence resulting in a steeper spectrum. As we 
already discussed in Section (4.2.1), the very flat spectra that we observe at the lowest flux 
states could be the result of the fact that in these cases the reflection component is prominent 
in the energy spectrum (affecting mainly the [7 — 10] keV band) so that the estimated slope 
determines mainly the shape of this component instead of the primary continuum. 

Apart from the average trend, the HRl—Hux plots also show signiflcant scatter around 
the best fltting lines. Therefore, spectral slope variations happen at constant flux as well. 
Since different slope values result in the same soft band ffux, these variations should be 
associated with luminosity changes. Apart from variability induced by the Compton-cooling 
process, other processes that affect the corona properties (like and the corona heating process 
for example) should also be variable. 
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6.2. Iron line and reflection component variability 

In NGC 5506 and MCG -6-30-15, the line's EW and the reflection parameter R vary 
proportional to each other (i.e. in the same way with the source flux). This implies that the 
material which emits the iron line is also responsible for the reflection hump. In the other 
two sources, R remains constant while EW decreases with the source flux. If this is the 
case, the line emitting and reflecting material should not be the same. However, it would 
be difficult to explain why matter that produces the iron line cannot produce the reflection 
hump as well (and vice versa). 

Changes in the line's EW and the reflection parameter it! could be the result of vari- 
ations in the geometry of the system. For example, if the sohd angle subtended by the 
reflector decreases, both the line's EW and the reflection parameter should also decrease. A 
relativistic outflow of the X-ray source could produce such an effect. However, at the same 
time, it should also be able to produce a softening of the primary spectrum and an increase 
in the total X-ray luminosity in order to explain the simultaneous increase of the soft band 
flux as both EW and R decrease. 

The observed EW and R variability behavior of the sources can be explained in a 
consistent way if the reprocessing material does not respond instantly to the continuum 
spectral variations. As the solid lines in Figures 4 and 5 show, the assumption of constant 
flux for the line and the reflection component during the spectral variations can account 
for the mean trends in the HR2 and HR3 color-flux diagrams. As the source luminosity 
and hence the continuum normalization increases, the soft band flux increases (see Figure 7) 
while both EW and R decrease because the flux of the line and reflection component remain 
constant. Since the average spectral properties of the sources are not identical (see Table 
3), the constant flux model curves can explain quite well the spectral variations although 
the best fitting lines in the color-fiux plots are not the same for all sources. Even in the 
cases where R remains constant during the fiux variations (e.g. NGC 5548) the model curves 
can still fit well the color-fiux diagrams. The reason is that the constant fiux and constant 
R model curves are quite similar; we need data with smaller uncertainty to discriminate 
between the two cases. 

A natural explanation for the fact that the reprocessor does not respond to the fast 
intrinsic variations is that it is located away from the central source, i.e. it is associated 
with the obscuring torus in the objects. Interestingly, this possibility, at least in the case 
of NGC 5548, is supported by a recent CHANDRA observation of this source which has 
revealed the existence of a narrow Fe K emission with a center energy of ~ 6.40 keV and 
EW ~ 130 eV (Yaqoob et al. 2001), similar to the average line's EW we find for this 
object (Table 3). As these authors comment, the fine should be produced in material that is 
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located away from the central source, in agreement with our results from the hardness ratios 
variations. 

Our results do not exclude reprocessing from material close to the central source (i.e. 
from the accretion disk). The HR2 color- flux diagram shows significant scatter around the 
best fitting lines. These variations could be the result of the slope variations that cause the 
scatter in the HRl color- flux plot (the [7 — 10] kcV band count rate contributes to both 
the HRl and HR2 ratios). However, the HR2 variations at constant flux could also be the 
result of fast variations of the iron line flux. 

Finally, we find a correlation between the spectral parameters of the average spectra of 
the individual sources in the sense that softer spectra show a lager reflection component R 
and line's EW value. This is similar to the correlation between R and V of Zdziarski et al. 
(1999) for different sources. However, we do not flnd the same correlation between R and F 
during the spectral variations within each source. Our results show clearly that the reflection 
parameter R is not increasing with increasing source flux (i.e. as the spectrum softens). The 
i?-flux correlation that we flnd is opposite to what we would expect if R was decreasing with 
decreasing F. Therefore, we do not conflrm the Zdziarski et al. (1999) correlation between 
V — R during the spectral variations within each source. If the bulk of reprocessing originates 
from neutral material located away from the central source and the T — R relationship is 
valid for different sources, then its origin is not obvious. It probably depends on a physical 
parameter that varies from one object to the other, like black hole mass or luminosity (which 
depends on the accretion rate as well). For example, in the more luminous sources the ffux of 
the soft photons could be stronger resulting in a stronger cooling of the plasma hence a softer 
spectrum. At the same time, the effective solid angle subtended by the reflector should also 
be larger, either because the reflector is located closer to the central source (which seems 
unlikely), or because the reflector is having a larger size. 



7. Conclusions 

We have used monitoring RXTE light curves to study the spectral variability in four 
Seyfert galaxies, namely NGC 4051, MCG -6-30-15, NGC 5506 and NGC 5548. The large 
number of observations over a long period (over ~ 3 years) allowed us to study their overall 
spectral variability behavior over many, different flux states. This is not possible to achieve 
with the normal observations of AGN which last typically less than a few days. 

The exposure time of each observation is not long enough to perform model fltting to 
the [3 — 15] keV band energy spectrum of the sources, but adequate to give us accurate 



-19- 



measurements of the flux in different energy bands. For tfiat reason, we computed fiardness 
ratios, we produced color-flux diagrams and we used the commonly used "power law plus 
reflection plus iron line" model to interpret our results. 

All sources show similar spectral variations. The variations in the HRl color-flux plot 
imply primary slope changes which affect significantly the soft band flux. The mean trends 
in the HR2 and HR3 color plots can be explained if wc assume that the flux of the line 

and of the reflection component remains roughly constant during the underlying continuum 
variations. A constant line and reflection component flux is naturally expected if the repro- 
cessing material is located away from the central source. The fact that the model under the 
assumption of constant line and reflection flux describes well the mean trend in the color- 
flux of all sources suggests that they operate in the same way. What are different is the 
amphtude of the slope variations (NGC 4051 shows the largest amphtude slope variability) 
and the average spectrum of each source. Interestingly, we find that sources with softer 
energy spectra show larger EW and R values as well. However, within each source, we find 
the opposite behavior. The reflection parameter R and the line's EW decrease (instead of 
increasing) with increasing source flux (and hence softening of the energy spectrum). There- 
fore, we conflrm the R — T relationship of Zdziarski et al. (1999) for the average spectra of 
the sources, but not for the spectral variations within each source. 

The signiflcant scatter around the best fltting lines in the HRl and HR2 color-flux 
diagrams shows that the spectral variability in AGN is more complicated than the mean 
trends of the color-flux plots suggest. The low amplitude slope variations at constant flux 
imply that Compton-cooling is not the only process that determines the properties of the 
active corona. The HR2 variations around the mean trend could imply low amplitude, line 
flux variations. If true, they could be the result of reprocessing from material near the central 
source. 

We thank the RXTE team for their operation of the satellite. lEP thanks the Osserva- 
torio Astronomico di Brera for hospitality. We also thank an anonymous referee for making 
helpful suggestions. Part of this work was done in the TMR research network 'Accretion 
onto black holes, compact stars and protostars' funded by the European Commission under 
contract number ERBFMRX-CT98-0195. 

APPENDIX 

FALSE CORRELATIONS IN THE COLOR-COLOR PLOTS 

As stated in the main text, the use of the same count rate in the determination of 
two colors can introduce misleading correlations in the color-color diagrams. To illustrate 
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this, in the top plot in Figure 8 we have plotted the HR2 versus HRl colors in the case of 
MCG -6-30-15. There appears to exist a strong correlation between the two ratios in the 
sense that HR2 increases with increasing HRl. The solid line in this plot shows the best 
fitting line to the color-color diagram (we find a best fitting slope of 1.24 ± 0.09). This is 
not consistent with the HR2 vs C (^-^keV)norm color-flux diagram of the same source which 
showed that HR2 remains constant, independent of the source flux state. Since the spectral 
slope is correlated with the source flux, we did not expect to see a correlation between HR2 
and HRl either. 

We performed a numerical experiment to verify whether the observed correlation be- 
tween HR2 and HRl is an artifact, caused by the inter-dependence of the two variables. 
First, we assumed 124 points (equal to the number of points in the MCG -6-30-15 hght curve) 
evenly distributed in the interval between 0.4 and 1.8 (i.e. the limits of the C'(3-5feey)„oT.„, val- 
ues of MCG -6-30-15) . Using the best fltting power law model line of the HRl vs C{z-bkeV)norm 
color-flux diagram, we computed the HRl values that correspond to these values, and hence 
the respective [7 — 10] keV band count rate. To compute the C^-jkeV keV values, we used 
the Cj-iokeV values and assumed a constant HR2 value (equal to the mean HR2 value listed 
in Table 1) at each flux state. Then we added appropriate random Gaussian "errors" to the 
simulated Cs^^keV, C's-TfeeV and Cy-iofceV points and computed the HRl and HR2 colors 
using the randomized simulated count rates. 

The simulated HRl and HR2 vs. C(^3^5keV)norrn color-flux diagrams arc shown in Figure 
8 (second and third plot from top). They arc similar to the respective diagrams in Figures 
3 and 4. The bottom plot in Figure 8 shows the HR2 vs HRl diagram using the simulated 
colors. Despite the fact that we had computed the C^-neV simulated values in such a way 
so that HR2 would remain constant at the same value irrespective of the source flux or 
spectral state, the simulated color-color plot shows a strong correlation between HR2 and 
HRl, as if larger values of HR2 correspond to larger values of HRl as well. In fact, the 
slope of the best fitting line to the simulated color-color plot (shown with a solid line in 
Figure 8) is 1.49 ± 0.17. The difference between this value and the best fitting slope value 
of the observed color-color diagram is 0.25 ± 0.17. This result shows that purely statistical 
effects are highly likely to reproduce the observed correlation between HRl and HR2. Since 
both ratios use C^^iokev in the numerator, when Cf^iokeV is higher than average due to a 
statistical fluctuation, then both HRl and HR2 will be higher on average, resulting in an 
almost one-to-one correlation between them. 
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Fig. 1.— The RXTE 2-10 keV monitoring light curves of NGC 4051, MCG -6-30-15, 

NGC 5506 and NGC 5548. The hght curves are normahzed to their weighted mean. The 
time axis is on days since the first observation of each source. Note the different samphng 
scheme throughout the observing period. Errors are included in the plotted points but are 
too small to be seen. 

Fig. 2. — The hardness ratio hght curves for the four sources. The sohd lines show the mean 
value of each light curve. Time is measured (in days) since first observation of each source. 
Note the different y— axis scaling in the case of NGC 4051. 

Fig. 3. — The HRl values plotted as a function of the normahzed [3 — 5] + [7 — 10] keV source 
count rate. The dashed hues show the best fitting power law to the data. Filled squares 
show the slope of a power law plus reflection model that corresponds to the respective [flux, 
HRl] values. These slope values were computed as described in Section 4.1 

Fig. 4. — The HR2 values plotted as a function of the normahzed [3 — 5] keV source count 
rate. Filled squares show the average HR2 values in various flux bins. In all plots, the 
dashed, dotted and dashed-dotted hues show the color-flux model curves in the case of a 
power law plus reflection plus iron hue with constant EW — 0, 250 and 500 eV respectively. 
Solid lines show the model color-flux curves in the case of a constant flux line. The model 
curves were computed as explained in Section 4.1 

Fig. 5. — The HR3 values plotted as a function of the normahzed [3 — 5] keV source count 
rate. Filled squares show the average HR3 values in various flux bins. In all plots, the 
dashed, dotted and dashed-dotted hues show the color-flux model curves in the case a power 
law plus reflection with constant R = 2,1 and respectively. Solid lines show the model 
color-flux curves in the case of a constant flux reflection. The model curves were computed 
as explained in Section 4.1 

Fig. 6. — The HRl values plotted as a function of the normalized [3 — 5] -|- [7 — 10] keV 
source count rate. As in Figure 3, the dashed lines show the best fltting power law to the 
data. Filled squares show the slope of the power law that corresponds to the respective 

[flux, HRl] values in the case when there exists a reflection component whose flux remains 
constant. These values were computed as described in Section 4.2.2. Note that, in this 
Figure, we have included the data from the low state of NGC 4051. 

Fig. 7. — Plot of the model, [1 — 300] keV flux as a function of the normahzed, soft band 
count rate. The peak-to-peak variation in the source luminosity is less than ~ 2 in all cases. 



Fig. 8.— Plot of the HR2 vs HRl color-color diagram in the case of MCG -6-30-15 (a). A 
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strong correlation is observed. The solid line shows the best fitting line to the data. The 
panels (b) and (c) show a plot of the HRl and HR2 vs normalized [3 — 5] keV count rate 
color-fiux diagrams for a set of simulated light curves (see Appendix). Panel (d) shows the 
HR2 vs HRl color-color diagram for the same set of simulated light curves. Although the 
light curves were constructed in such a way so that HR2 is constant, there exists a strong 
correlation between HR2 and HRl (similar to the correlation between the observed HRl 
and HR2 values of MCG -6-30-15 shown in panel (a)). This correlation is caused by random, 
statistical effects. 
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Namc 



H Rl (Trrns H R1 O'rms H R3 (T„ 



NGC 4051 0.37 0.17 0.57 0.06 0.58 0.02 

MCG -6-30-15 0.41 0.06 0.61 0.04 0.60 0.05 

NGC 5506 0.51 0.03 0.64 0.03 0.62 0.04 

NGC 5548 0.47 0.04 0.67 0.02 0.61 



Table 1: The mean and (Trms values of the HRl, HR2, and HR3 light curves shown in 
Figure 2. The error on the mean HR values is less than 0.01. 



Object Name 


HRl/ C(3-5keV+7-10keV)„orm 

bi (xVdof) 


hi (xVdof) 


HR2,/ C(3_5fcey)„o,^ 
hi (xVdof) 


NGC 4051 
MCG -6-30-15 
NGC 5506 
NGC 5548 


-0.30 ±0.02(284/123) 
-0.16 ±0.02(256/122) 
-0.10 ±0.01(293/125) 
-0.13 ±0.01(290/132) 


-0.05 ±0.01(189/123) 
+0.02 ±0.02(190/122) 
-0.01 ±0.01(193/125) 
+0.01 ±0.01(191/132) 


-0.17 + 0.02(124/123) 
-0.13 ±0.03(132/122) 
-0.10 ±0.01(130/125) 
-0.06 ±0.02(110/132) 



Table 2: Best power law model fitting results to the color- fiux plots shown in Figures 3, 4 
and 5. 
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Name 


(r) 


{R) 


(EW) (eV) 


NGC 4051 


2.06 


1.7 


310 


MCG -6-30-15 


2.05 


1.2 


180 


NGC 5506 


2.11 


1.4 


230 


NGC 5548 


1.88 


0.8 


130 



Table 3: The mean spectral index, reflection normalization R and iron line equivalent width 
(EW) of the sources in our sample, computed using the mean hardness ratio values and a 
power law plus reflection and a narrow Gaussian line spectral model (see Section 4.2.1). 



Name HR2/ C(^3-5keV)„orn, HR3/ C(^3_5keV). 





bcEW 


bcLF 


bcR 


bcRF 


NGC 4051 


-0.14 


-0.03 


-0.13 


-0.15 


MCG -6-30-15 


-0.08 


-0.02 


-0.08 


-0.11 


NGC 5506 


-0.04 


-0.02 


-0.05 


-0.09 


NGC 5548 


-0.06 


-0.03 


-0.06 


-0.09 



Table 4: The slope of the model curves shown in Figures 4 and 5. bcEW and bcLF represent the 
slope of the model curves under the assumption of constant line's EW and flux respectively. 
Similarly, hcR and hcRF denote the slope of the model curves in the case of constant it! and 
constant reflection flux respectively. 



